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We study the anisotropic boundary conditions for the dilute 0{n) loop model with the methods of 2D 
quantum gravity. We solve the problem exactly on a dynamical lattice using the correspondence with 
a large N matrix model. We formulate the disk two-point functions with ordinary and anisotropic 
boundary conditions as loop correlators in the matrix model. We derive the loop equations for these 
correlators and find their explicit solution in the scaling limit. Our solution reproduces the boundary 
phase diagram and the boundary critical exponents obtained recently by Dubail, Jacobsen and Saleur, 
except for the cusp at the isotropic special transition point. Moreover, our solution describes the bulk 
and the boundary deformations away from the anisotropic special transitions. In particular it shows 
how the anisotropic special boundary conditions are deformed by the bulk thermal flow towards the 
dense phase. 
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1 Introduction 



The boundary critical phenomena appear in a large spectrum of disciplines of the contemporary the- 
oretical physics, from solid state physics to string theory. The most interesting situation is when 
the boundary degrees of freedom enjoy a smaller symmetry than those in the bulk. In this case one 
speaks of surface anisotropy. The D-branes in string theory are perhaps the most studied example 
of such anisotropic surface behavior. Another example is provided by the ferromagnets with surface 
exchange anisotropy, which can lead to a critical and multi-critical anisotropic surface transitions. An 
interesting but difficult task is to study the interplay of surface and the bulk transitions and the related 
multi-scaling regimes. 

For uniaxial ferromagnets as the Ising model, there are four different classes of surface transi- 
tions: the ordinary, extraordinary, surface, and special transitions [1]. This classification makes sense 
also for spin systems with continuous 0{n) symmetry. It was predicted by Diehl and Eisenriegler 
[2], using the e-expansion and renormalization group methods, that the effects of surface anisotropy 
can be relevant near the special transitions of the d-dimensional 0(n) model. These effects lead to 
'anisotropic special transitions' with different critical exponents. 

Recently, an exact solution of the problem for the 2-dimensional 0(n) model was presented by 
Dubail, Jacobsen and Saleur [3] using its formulation as a loop model [4, 5]. Using an elaborate 
mixture of Coulomb gas, algebraic and Thermodynamic Bethe Ansatz techniques, the authors of [3] 
confirmed for the dilute 0{n) model the phase diagram suggested in [2] and determined the exact 
scaling exponents of the boundary operators. A review of the results obtained in [3], which is acces- 
sible for wider audience, can be found in [6]. Their works extended the techniques developped by 
Jacobsen and Saleur [7, 8] for the dense phase of the 0(n) loop model. 

In this paper we examine the bulk and the boundary deformations away from the anisotropic 
special transitions in the two-dimensional 0(n) model. In particular, we address the question how the 
anisotropic boundary transitions are influenced by the bulk deformation which relates the dilute and 
the dense phases of the O(n) model. To make the problem solvable, we put the model on a dynamical 
lattice. This procedure is sometimes called 'coupling to 2D gravity' [9, 10]. The sum over lattices 
erases the dependence of the correlation functions on the coordinates, so they become 'correlation 
numbers'. Yet the statistical model coupled to gravity contains all the essential information about the 
critical behavior of the original model such as the qualitative phase diagram and the conformal weights 
of the scaling operators. When the model is coupled to gravity, the bulk and boundary flows, originally 
driven by relevant operators, becomes marginal, the Liouville dressing completing the conformal 
weights to one. This necessitates a different interpretation of the flows. The UV and the IR limits 
are explored by taking respectively large and small values of the bulk and boundary cosmological 
constants. Our method of solution is based on the mapping to the 0(n) matrix model [11, 12] and 
on the techniques developed in [13, 14]. Using the Ward identities of the matrix model, we were able 
to evaluate the two-point functions of the boundary changing operators for finite bulk and boundary 
deformations away from the anisotropic special transitions. 

The paper is structured as follows. In Sect. 2 we summarize the known results about the boundary 
transitions in the 0(n) loop model. In Sect. 3 we write down the partition function of the boundary 
0(n) model on a dynamical lattice. In particular, we give a microscopic definition of the anisotropic 
boundary conditions on an arbitrary planar graph. In Sect. 4 we reformulate the problem in terms of 
the O(n) matrix model. We construct the matrix model loop observables that correspond to the disk 
two-point functions with ordinary and anisotropic special boundary conditions. In Sect. 5 we write a 
set of Ward identities (loop equations) for these loop observables, leaving the derivation to Appendix 
A. We are eventually interested in the scaling limit, where the volumes of the bulk and the boundary 
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of the planar graph diverge. This limit corresponds to tuning the bulk and the boundary cosmological 
constants to their critical values. In Sect. 6 we write the loop equations in the scaling limit in the form 
of functional equations. From these functional equations we extract all the information about the 
bulk and the boundary flows. In particular, we obtain the phase diagram for the boundary transitions, 
which is qualitatively the same as the one suggested in [3, 6], apart of the fact that we do not observe 
a cusp near the special point. In Sect. 7 we derive the conformal weights of the boundary changing 
operators. All our results concerning the critical exponents coincide with those of [3, 6]. In Sect. 8 we 
find the explicit solution of the loop equations in the limit of infinitely large planar graph. The solution 
represents a scaling function of the coupling for the bulk and boundary perturbations. The endpoints 
of the bulk and the boundary flows can be found by taking different limits of this general solution. 
The boundary flows relate the anisotropic special transition with the ordinary or with the extraordinary 
transition, depending on the sign of the perturbation. The bulk flow relates an anisotropic boundary 
condition in the dilute phase with another anisotropic boundary condition in the dense phase. For the 
rational values of the central charge, the boundary conditions associated with the endpoints of the bulk 
flow match with those predicted in the recent paper [15] using perturbative RG techniques. 



2 The boundary 0(n) model on a regular lattice: a summary 

The 0(n) model [4] is one of the most studied statistical models. For the definition of the partition 
function see Sect. 3. The model has an equivalent description in terms of a gas of self- and mutually 
avoiding loops with fugacity n. The partition function of the loop gas depends on the temperature 
coupling T which controls the length of the loops: 

Zq(, ) = n #[ lo °P s ] (]yy)#[ links occupied by loops] ^ 1) 

loops 

In the loop gas formulation of the 0(n) model, the number of flavors n can be given any real value. 
The model has a continuum transition if the number of flavors is in the interval [—2, 2]. Depending 
on the temperature coupling T the model has two critical phases, the dense and the dilute phases. At 
large T the loops are small and the model has no long range correlations. The dilute phase is achieved 
at the critical temperature [16, 17] 



T c = y 2 + V2 - n (2.2) 

for which the length of the loops diverges. If we adopt for the number of flavors the standard 
parametrization 

n = 2cos(tt6>), < 6 < 1, (2.3) 

then at the critical bulk temperature the 0(n) model is described by (in general) a non-rational CFT 
with central charge 

e 2 

^ = l~G—. (2.4) 

The primary operators <J> r s in such a non-rational CFT can be classified according to the generalized 
Kac table for the conformal weights, 

(r,-.)' -b-l)» 
4ff 
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Unlike the rational CFTs, here the numbers r and s can take non-integer values. 

When T < T c , the loops condense and fill almost all space. This critical phase is known as the 
dense phase of the loop gas. The dense phase of the 0(n) is described by a CFT with lower value of 
the central charge, 
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^ = 1-6^- (2-6) 
The generalized Kac table for the dense phase is 

(r~g-s) 2 -(g-l) 2 ,1 
h rs = ^ , 9 = — Q . (2-7) 

Most of the exact results for the dense and the dilute phases of the 0(n) model were obtained by 
mapping to Coulomb gas [5, 16]. 

The boundary 0(n) model was originally studied for the so called ordinary boundary condition, 
where the loops avoid the boundary as they avoid themselves. The ordinary boundary condition is 
also referred as Neumann boundary condition because the measure for the boundary spins is free. The 
boundary scaling dimensions of the L-leg operators Sl, realized as sources of L open lines, were 
conjectured for the ordinary boundary condition in [18] and then derived in [19]: 

(Ord\S L \Ord) ->$f +L .! (dilute phase) 

{Ord\S L \Ord) -^<S>f <1+L (dense phase). (2.8) 

Another obvious boundary condition is the fixed, or Dirichlet, boundary condition, which allows, 
besides the closed loops, open lines that end at the boundary [20]. 2 The dimensions of the L-leg 
boundary operators with Dirichlet and Neumann boundary conditions were computed in [21, 22] by 
coupling the model to 2D gravity and then using the KPZ scaling relation [23, 24]: 

(Ord\S L \Dir) ->• $f/ 2+L ,o (dilute phase) 

{Ord\S L \Dir) -» $J 1/2+L (dense phase). (2.9) 

From the perspective of the boundary CFT, these operators are obtained as the result of the fusion 
of the L-leg operator and a boundary-condition-changing (BCC) operator, introduced in [25], which 
transforms the ordinary into Dirichlet boundary condition. 

Recently it was discovered that the 0(n) loop model can exhibit unexpectedly rich boundary crit- 
ical behavior. Jacobsen and Saleur [7, 8] constructed a continuum of conformal boundary conditions 
for the dense phase of the 0(n) loop model. The Jacobsen-Saleur boundary condition, which we 
denote shortly by JS, prescribes that the loops that touch the boundary at least once are taken with 
fugacity y / n, while the loops that do not touch the boundary are counted with fugacity n. The 
boundary parameter y can take any real value. The JS boundary conditions contain as particular cases 
the ordinary (y = n) and fixed (y = 1) boundary conditions for the 0(n) spins. Jacobsen and Saleur 
conjectured the spectrum and the conformal dimensions of the L-leg boundary operators separating 
the ordinary and the JS boundary conditions. These conformal dimensions were subsequently veri- 
fied on the model coupled to 2D gravity in [13, 14]. If y is parametrized in the 'physical' interval 
< y < n as 

sin7r(r + 1)9 , .„ 

y = , V V (l<r< 1/0-1), (2.10) 

sin irrd 



- In the papers [20, 21, 22] the loop gas was considered in the context of the SOS model, for which the Dirichlet and 
Neumann boundary conditions have the opposite meaning. 
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then the BCC operator transforming the ordinary into JS boundary condition is identified as the diag- 
onal operator <£ r r . Note that r here is not necessarily integer or even rational. The L-leg operators 
with L > 1 fall into two types. The operator creates open lines which avoid the JS boundary. The 
operator creates open lines that can touch the JS boundary without restriction. The two types of 
L-leg boundary operators are identified as 

{Ord\Sf\JS) -^$r,r±L (dense phase). (2.11) 

The general case of two different JS boundary conditions with boundary parameters y\ and j/2 was 
considered for regular and dynamical lattices respectively in [26] and [27]. 

The JS boundary condition was subsequently adapted to the dilute phase by Dubail, Jacobsen and 
Saleur [3]. The authors of [3] considered the loop gas analog of the anisotropic boundary interaction 
studied previously by Diehl and Eisenriegler [2], which breaks the symmetry as 

0(n) -> 0(n«) x 0(n< 2 )), n« + n( 2 > = n. (2.12) 

The boundary interaction depends on two coupling constants, A(i) and A(2), associated with the two 
unbroken subgroups. In terms of the loop gas, the anisotropic boundary interaction is defined by 
introducing loops of two colors, (1) and (2), having fugacities respectively nW and n^. Each time 
when a loop of color (a) touches the JS boundary, it acquires an extra factor A(«) . We will call this 
boundary condition dilute Jacobsen-Saleur boundary condition, or shortly DJS , after the authors of 
[3]. 

Let us summarize the qualitative picture of the surface critical behavior proposed in [3]. Consider 
first the isotropic direction A(i) = A(2) = A. In the dilute phase one distinguishes three different 
kinds of critical surface behavior: ordinary, extraordinary and special. When A = 0, the loops in 
the bulk almost never touch the boundary. This is the ordinary boundary condition. When w — > oo, 
the most probable loop configurations are those with one loop adsorbed along the boundary, which 
prevents the other loops in the bulk to touch the boundary. The adsorbed loop plays the role of a 
boundary with ordinary boundary condition. This is the extraordinary transition. In terms of the 0(n) 
spins, ordinary and the extraordinary boundary conditions describe respectively disordered and or- 
dered boundary spins. 3 The ordinary and the extraordinary boundary conditions describe the same 
continuous theory except for a reshuffling of the boundary operators. The L-leg operator with or- 
dinary/extraordinary boundary conditions will look, when L > 1, as the (L — l)-leg operator with 
ordinary/ordinary boundary conditions, because its rightmost leg will be adsorbed by the boundary. 
The 0-leg operator with ordinary/extraordinary boundary condition will look like the l-leg operator 
with ordinary/ordinary boundary condition because one of the vacuum loops will be partially adsorbed 
by the extraordinary boundary and the part which is not adsorbed will look as an open line connecting 
the endpoints of the extraordinary boundary. 

The ordinary and the extraordinary boundary conditions are separated by a special transition, 
which happens at some A = A c and describes a conformal boundary condition. For the honeycomb 
lattice the special point is known [28] to be at 

A c = (2 - n)~ 1/2 T c 2 (honeycomb lattice). (2.13) 

At the special point the loops touch the boundary without being completely adsorbed by it. The special 
transition exists only in the dilute phase, because in the dense phase the loops already almost surely 

3 The spontaneous ordering on the boundary does not contradict the Mermin- Wagner theorem. In the interval 1 < n < 2 
the target space, the (n — 1) -dimensional sphere, has negative curvature and thus resembles a non-compact space. 
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touch the boundary. The only effect of having small or large A is the reshuffling of the spectrum of 
L-leg operators, which happens in the same way as in the dilute phase. 

In the anisotropic case, A(i) 7^ A(2), there are again three possible transitions: ordinary, extraor- 
dinary and special. When A(i) and A(2) are small, we have the same ordinary boundary condition as 
in the isotropic case. In the opposite limit, where A(i) and A(2) are both large, the boundary spins be- 
come ordered in two different ways, depending on which of the two couplings prevails. If A(i) > A(2) , 
the (l)-type components order, while the (2)-type components remain desordered, and vice versa. 
The (A(i) , A(2) ) plane is thus divided into three domains, characterized by disordered, (l)-ordered and 
(2)-disordered, and (2)-ordered and (1) -disordered, which we denote respectively by Ord, Ext^ 
and Ext( 2 ) - The domains Ext(\\ and Ext( 2 ) are separated by the isotropic line starting at the special 
point and going to infinity. This is a line of first order transitions because crossing it switches from 
one ground state to the other. The remaining two critical lines, Ord/ Ext^ and Ord/ Ext^, are the 
lines of the two anisotropic special transitions, AS^ and AS( 2 y K was argued in [2, 3, 6], using 
scaling arguments, that the lines AS^ and AS( 2 ) join at the point Sp = (A c , A c ) in a cusp-like shape. 
The model was solved in [3] for a particular point on AS^y. 



1 - n( 2 ) + VI - nWnP) 
A(i) = H 



V2 
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A( 2) = 1 + . (2.14) 

a/2 — n 

In terms of the loop gas expansion, the anisotropic special transitions are obtained by critically 
enhancing the interaction with the boundary of the loops of color (1) or (2). The boundary CFT's 
describing the transitions ASn\ and AS( 2 ) were identified in [3]. A convenient parametrization of 
nW and on the real axis is 4 

sinf7r(r - 1)01 sin \n (r + 1)01 , 

n(» = L .\ -/ J , = l . \ V J (0 < r < 1/0) . (2.15) 

sm[7rrc/J sin[7rrc/J 

The loop model has a statistical meaning only if both fugacities are positive, which is the case when 
1 <r< 1/0 — 1. With the above parametrization, the BCC operators (AS^\Ord) and (AS (2) \Ord) 
are argued to be respectively <& pr and $f T+1 . More generally, one can consider the L-leg boundary 

operators, S^ and S^\ which create L open lines of color respectively (1) and (2). The Kac labels 
of these operators were determined in [3] as follows, 

(Ord\S?\AS {1) ) - $f +Lr , (Ord\S?\AS {1) ) *f_ i>p) 
(Ordl^l^a)) -$ r 3 +i , r+1 , (Ord\Sf\AS {2) ) - *f_ £>P+1 . (2.16) 

In the vicinity of the special transitions the theory is argued to be described by a perturbation of the 
boundary CFT by the boundary operator $>f 3 in the isotropic direction and by in the anisotropic 
direction. 



4 The correspondence between our notations and the notations used in [6] is {nW , n (2 > , r, A(i) , A<2) }here 



{n2,m, r, W2, Wl}there 
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3 The boundary 0(n) model on a dynamical lattice 



3.1 Anisotropic boundary conditions for the O(n) model on a planar graph 

The two-dimensional 0(n) loop model, originally denned on the honeycomb lattice [4], can be also 
considered on a honeycomb lattice with defects, such as the one shown in Fig. la. The lattice repre- 
sents a trivalent planar graph T. We define the boundary dT of the graph by adding a set of extra lines 
(the single lines in the figure) which turn the original planar graph into a two-dimensional cellular 
complex. The local fluctuating variable is an 0(n) classical spin, that is an ?7,-component vector S(r) 
with unit norm, associated with each vertex r G T, including the vertices on the boundary dT. The 
partition function of the 0(n) model on the graph T depends on the coupling T, called temperature, 
and is defined as an integral over all classical spins, 

Z 0(n) (T;T)= I \{[dS{r)]\{ (l + A S(r) ■ S(r')) , (3.1) 

rer (rr 1 ) 

where the product runs over the lines {rr') of the graph, excluding the lines along the boundary. The 
0(n) -invariant measure [dS] is normalized so that 

J [dS] S a S b = 5 a>b . (3.2) 

The partition function (3.1) corresponds to the ordinary boundary condition, in which there is no 
interaction along the boundary. 

Expanding the integrand as a sum of monomials, the partition function can be written as a sum 
over all configurations of self-avoiding, mutually-avoiding loops as the one shown in Fig. lb, each 
counted with a factor of n: 

Z 0{n) {T-T)= ^ T' len ^ th n #l0 °P s . (3.3) 

loops on r 

The temperature coupling T controls the length of the loops. The advantage of the loop gas represen- 
tation (3.3) is that it makes sense also for non-integer n. In terms of loop gas, the ordinary boundary 
condition, which we will denote by Ord, means that the loops in the bulk avoid the boundary as they 
avoid the other loops and themselves. 

The Dirichlet boundary condition, was originally defined for the dense phase of the loop gas 
[20, 21, 22] and requires that an open line starts at each point on the boundary. The dilute version 
of this boundary condition depends on an adjustable parameter, which controls the number of the 
open lines. In terms of the 0{n) spins the Dirichlet boundary condition is obtained by switching on 
a constant magnetic field B acting on the boundary spins. This modifies the integration measure in 
(3.1) by a factor 

JJ (l + B ■ S{r)) , (3.4) 
rear 

where the product goes over all boundary sites r. The loop expansion with this boundary measure 
contains open lines having both ends at the boundary, each weighted with a factor B 2 . 

The dilute anisotropic (DJS) boundary condition is defined as follows. The n components of the 
0(n) spin are split into two sets, (1) and (2), containing respectively nW and components, with 
n (i) _|_ n (2) _ n This leads to a decomposition of the 0(n) spin as 

S = S(i) + S(2), • 5(2) = 0. (3.5) 
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Figure 1 : a) A trivalent planar graph T with a boundary b) A loop configuration on T for the ordinary 
boundary condition. The loops avoid the boundary as they avoid themselves 




Figure 2: A loop configuration for the JS boundary condition. The loops in the bulk (in red) have fugacity 
n, while the loops that touch the boundary (in blue and green) have fugacities nw or n< 2 > depending on their 
color. 

The DJS boundary condition is introduced by an extra factor associated with the boundary links, 

H (l + Awi)(r)'i)(r')). (3.6) 

(rr')Gdr a=l,2 

This boundary interaction is invariant under the subgroup of independent rotations of Sn\ and S^ 2 )- 
The boundary term changes the loop expansion. The loops are now allowed to pass along the boundary 
links as shown in Fig. 2. We have to introduce loops of two colors, (1) and (2), having fugacities 
respectively raM and . A loop of color (a) that visits J\f boundary links acquires an additional 
weight factor . For the loops that do not touch the boundary, the contributions of the two colors 
sum up to + n' 2 > = n and we obtain the same weight as with the ordinary boundary condition. 

3.2 Coupling to 2D discrete gravity 

The disk partition function of the 0{n) model on a dynamical lattice is defined as the expectation value 
of (3.1) in the ensemble of all trivalent planar graphs T with the topology of the disk. The measure 
depends on two more couplings, fl and fig, called respectively bulk and boundary cosmological con- 
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stants, 5 associated with the volume |T| = #(cells) and the boundary length \dT\ = ^(external lines). 
The partition function of the disk is a function of p, and p, b and is denned by 




|sr| 

Z 0( JT;T). (3.7) 



3.3 Two-point functions of the L-leg boundary operators 

Our aim is to evaluate the boundary two-point function of two L-leg operators separating ordinary 
and anisotropic boundary conditions. The L-leg operator Sl is obtained by fusing L spins with flavor 
indices a±, . . . , a n € {1, . . . , n}. In terms of the loop gas, the operator Sl creates L self and mutually 
avoiding open lines. We would like to exclude configurations where some of the lines contract among 
themselves. This can be achieved by taking the antisymmetrized product 

5 L ~ detSo.fo), (3.8) 

where n, . . . , rj_, are L consecutive boundary vertices of the planar graph T, and we put the label L 
instead of writing its dependence on oi, . . . , cll, explicitly. The two-point function of the operator Sl 
is evaluated as the partition function of the loop gas in presence of L open lines connecting the points 
{ri} and {r',i}. The open lines are self- and mutually avoiding, and are not allowed to intersect the 
vacuum loops. 

Since the DJS boundary condition breaks the 0(n) symmetry into 0(nW) x 0(n^), there are 
two inequivalent correlation functions of the L-leg operators with Or d/ DJS boundary conditions. 
Indeed, the insertion of S a has different effects depending on whether a belongs to the O(nW) or the 
0(n( 2 ') sectors. In the first case the open line created by S a acquires a factor A(i) each time it visites 
a boundary link. In the second case, the factor is A(a). Therefore the boundary spin operators (3.8) 
with Ord/DJS boundary conditions split into two classes, 

S[ 1} ~ det [S^r,)] , ai,...,a L G (1), 

4 2) ~ det[S ai ( rj )], oi,...,ol€(2). (3.9) 
We denote the corresponding boundary two-point functions respectively by and D L 2 ' . 

4 Mapping to the O(n) matrix model 

The 0(n) matrix model [11, 12] generates planar graphs covered by loops in the same way as the 
one-matrix models considered in the classical paper [29] generate empty planar graphs. The model 
involves the hermitian N x TV matrices M and Y a , where the flavor index a takes n values. The 
partition function is given by an 0(n)-invariant matrix integral 

Z N (T) ~ J dMd n Ye-^ M2+ ^ 2 ^ M3 ~ M ^). (4.1) 

This integral can be considered as the partition function of a zero-dimensional QFT with Feynman 
rules given in Fig. 3, where we used the 't Hooft double-line notations. The graphs made of such 

5 We used bars to distinguish from the bulk and boundary cosmological constants in the continuum limit. 
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a i 



1//3 1/0T 15 (3 

Figure 3: Feynman rules for the 0(n) matrix model 



i J 



ik kJ 



' — II — / 



ik kj 



^I^b A ( i)///| A( 2 )///| 

Figure 4: The constituents of the DJS boundary generated by M, Y?js and Y? 2 %: a non-occupied site, a 
sequence of two sites visited by a loop of color (1), and a sequence of two sites visited by a loop of color (1) 



double-lined propagators are known as fat graphs. The 'vacuum energy' of the matrix model repre- 
sents a sum over connected fat graphs, which can be also considered as discretized two-dimensional 
surfaces of all possible genera. As the action is quadratic in the matrices Y a , their propagators arrange 
in closed loops carrying a flavor a. The sum of all Feynman graphs with given connectivity can be 
viewed as the sum over all configurations of self and mutually avoiding loops on a given discretized 
surface. The weight of each loop is given by the product of factors 1/T, one for each link, and the 
number of flavors n. We are interested in the large N limit 

N -c oo, j3/N = fl 2 (fixed) , (4.2) 

in which only fat graphs of genus zero survive [30]. 

The basic observable in the matrix model is the resolvent 

^-K"^)' <4 - 3) 

evaluated in the ensemble (4.1). The resolvent is the one-point function with ordinary boundary 
conditions and is related to the disk partition function by W = —dp B U. 

The one-point function with Dirichlet boundary condition is obtained by adding a term B- Y which 
expresses the coupling with the magnetic field on the boundary. This leads to a more-complicated 
resolvent 

R(fl B ,B 2 ) = l(tr L__.\. (4.4) 

In order to include the anisotropic boundary conditions in this scheme, we decompose the vector 
Y into a sum of an raW -component vector Y(i) and an n< 2 ) -component vector Y( 2 ) as in (3.5): 

Y = Y (1) +Y (2) , Y (1) -Y (2) =0. (4.5) 



9 



The one-point function with ordinary and DJS boundary conditions is given by the resolvent 

H(ji B , A ( i) , A (2) ) = \ ( tr — — ^ ) . (4.6) 

The two extra terms are the operators creating boundary links containing segments of lines of type (1) 
and (2), as shown in Fig.4. Each such operator created two boundary sites, hence the factor \/J1b- 

The matrix integral measure becomes singular at M = T/2. We perform a linear change of the 
variables 

M = T(i + X), (4.7) 

which sends this singular point to X = 0. After a suitable rescaling of Y and f5, the matrix model 
partition function takes the canonical form 

Z N ~ J d X d n Ye ^[-V(X) + XY 2 ] ( (4 8) 

where V(x) is a cubic potential 

3 

V{x) = £ ~ x3 = ~ | (* + l) 3 + \{* + h) 2 - ( 4 - 9 ) 



3=0 

We introduce the spectral parameter x which is related to the lattice boundary cosmological constant 

b y 

fL B =T(x + l). (4.10) 
Now the one-point function with ordinary boundary condition is 

W(x) = - (trW(x)) , (4.11) 

P 

where the matrix 

W0r) d ^'— (4.12) 
x — X. 

creates a boundary segment with open ends. In the following we will call x boundary cosmological 
constant. We also redefine the boundary couplings X( a ) in (4.6) as 

A( Q ) — > A( a ) /2b, a = 1,2. (4.13) 

Then the operator that creates a boundary segment with DJS boundary condition is 

H(y) = l — — . (4.14) 

y - X - A(i) Y (1) - A(2)Y (2) 

The boundary L-leg operators are represented by the antisymmetrized products 

def 

Sl = Y ai Y a2 • • • Y aL ± permutations . (4.15) 

The boundary two-point function of the L-leg operator with Ord/Ord boundary conditions is given 
by the expectation value 

D L ( Xl , x 2 ) = j 3 <tr[W(x 1 )S L W(a; 2 )S L ] > . (4.16) 
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The role of the operators W(cci) and W(x2) is to create the two boundary segments with boundary 
cosmological constants respectively x\ and x%. The two insertions Sl generate L open lines at the 
points separating the two segments. It is useful to extend this definition to the case L = 0, assuming 
that Sq is the boundary identity operator. In this simplest case the expectation value (4.16) is evaluated 
instantly as 

D (xi,x 2 ) = iMWOrxJWM) = W ^~ W ^\ (4.17) 

p X\ — X2 

The two-point functions (4.16) for L > 1 were computed in [21, 22]. 

In the case of a DJS boundary condition, the matrix model realization of the two types of bound- 
ary L-leg operators is given by the antisymmetrized products (4.15), with the restrictions on the com- 
ponents as in (3.9). The boundary two-point functions with Ord/DJS boundary conditions are 
evaluated by the expectation values 

D (x,y) = ^tr[W(s)H(iO]) (4.18) 

and for L > 1, 

D®(x,y) = i(tr[W( a ;)si 1) H(y)S? ) ]), (4.19) 

Df{x,y) = i(tr[W(x)sfH(y)si 2) ]). (4.20) 

Apart of the DJS boundary parameter nW and the boundary couplings A(i), A(2) for the second seg- 
ment, they depend on the boundary cosmological constants x and y associated with the two segments 
of the boundary. 



5 Loop equations 

Our goal is to evaluate the two-point functions (4.19) in the continuum limit, when the area and 
the boundary length of the disk are very large. They will be obtained as solution of a set of Ward 
identities, called loop equations, which follow from the translational invariance of the integration 
measure in (4.1), and in which the n' 1 ) enters as a parameter. The solutions of the loop equations are 
analytic functions of nW which can take any real value. We will restrict our analysis to the 'physical' 
case < n* 1 ) < n, when the correlation functions have good statistical limit. Here we summarize 
the loop equations which will be extensively studied in following sections. The proofs are given in 
Appendix A. 

5.1 Loop equation for the resolvent 

The loop equation for the resolvent is known [31], but we nevertheless recall it here in order to set up 
a self-contained description of the method. The resolvent W(x) splits into a singular part w(x) and a 
polynomial W Teg (x): 

Hpf 

W(x) = W reg (x) + w(x) . (5.1) 
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The regular part is given by 

w < ^ 2V'(x)-nV'(-x) 2 
W veg {x) = ^ = -a - ciiX - a 2 x 



do 



a i 



2V'(x) - 


- nV'(-x) 


4 


— n 2 


91 


T-2 


2 + n 


4(2 + n) 


92 


T-l 


2-n 


2-n 


93 


T 


2 + n 


2 + n' 



" 2 = -2T^ = - 

The function w(x) satisfies a quadratic identity 

w 2 (x) + w 2 (-x) + nw(x)w(-x) = A + Bx 2 + Cx A . (5.3) 



The coefficients A, B, C as functions of T, a and W\ = (trX) can be evaluated by substituting the 
large-x asymptotics 

w (x) = -W rcg (x) + ^— + ^ + 0(x~ 3 ) (5.4) 



in (5.3). The solution of the loop equation (5.3) with the asymptotics (5.4) is given by a meromorphic 
function with a single cut [a, b] on the first sheet, with a < b < 0. This equation can be solved by an 
elliptic parametrization and the solution is expressed in terms of Jacobi theta functions [32]. 

5.2 Loop equations for the boundary two-point functions with Ord/DJS boundary 
conditions 

The two-point correlators (4.16) with ordinary boundary conditions are known to satisfy the integral 
recurrence equations [21, 22] 

D L+1 = W*D L . (5.5) 

The "^-product" is defined for any pair of meromorphic functions, analytic in the right half plane 
Re(x) > and vanishing at infinity, 

r , u . def f dx' f(x) - f(x') , 

[/**](*) = -f^- x _ x , 9(-x), (5-6) 
C- 

with the contour C_ encircling the left half plane Rex < 0. These equations actually hold for a more 
general set of two-point correlators, which have ordinary boundary condition on one segment and an 
arbitrary boundary condition on the other segment [13]. Thus the boundary two-point functions (4.19) 
and (4.20) for L > 1 satisfy the same recurrence equations 

DM, = W*D%\ 

D {2) +1 = W*D {2) . (5.7) 

Using the recurrence relation, the correlation functions of the L-leg operators can be obtained recur- 
sively from those of the one-leg operators and 
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The correlator Dq, defined by (4.18), and the correlators and which we normalize as 
D[ a \x,y) = ^^(tr[W(x)YWH(j,)Yi Q )]) (a = 1,2), (5.8) 

can be determined by the following pair of bilinear functional equations, derived in Appendix A. In 
order to shorten the expressions, here and below we use the shorthand notation 

F{x) = F(-x) . (5.9) 

The two equations then read 



W -H + D (x-y) + nWA ( .) ((\( a) D -l)D[ a) + WDo) = 0, (5. 

o=l,2 ^ ' 



10) 



P + D (H-V' + W + nW) + J2 nW(\ {a) D -l)D[ a) = 0. (5.11) 

a=l,2 



The second equation involves an unknown linear function of x: 

P(x,y) =* 1 ^ V ' {X) X ~ x (X) h ^)) = 92H(y) + gsH^y) + xg 3 H(y). (5.12) 

Equations (5. 10) and (5.11) can be solved in favor of or If we define 

A™ d = f X W D -1, 

= (A ( i) - A( 2 ))n(i){A (1)J Di 1) +W} - A (2) (W -V 7 + H) -x-y, (5.13) 
= XwX^P + (A(i) + \ (2) )H -x + y- (A(i) - A (2) ) (W + nWW) - \(2)V , 

and similarly for A^ 2 \B^ 2 \C^ 2 \ with an obvious exchange (1) <-> (2), then (5.10) and (5.11) take 
the following factorized form: 



A {a) B( a ) = C {a) (a = 1,2). (5.14) 

Equations (5.14) are the main instrument of our analysis of the DJS boundary conditions. 
It is convenient to define the functions Dq by 



D { « ] = = f ° = (« = 1,2). (5.15) 

1 — A(a)L>Q 



In Appendix A we show that with this definition the recurrence equations (5.7) hold also for L = 0. 
The equation for L = is a consequence of (5.14). 



5.3 Loop equation for the two-point function with Ord/Dir boundary conditions 

On the flat lattice, the DJS boundary condition with = 1 is equivalent, for a special choice of 
the boundary parameters, to the Dirichlet boundary condition defined by the boundary factor (3.4). In 
order to make the comparison on the dynamical lattice, we will formulate and solve the loop equation 
for the two-point function of the BCC operator with ordinary/Dirichlet boundary conditions. 
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Assume that the magnetic field points at the direction a = 1. Then the correlator in question is 
given in the matrix model by the expectation value 

n(x, y) = \ ( trW(x)R(i/) ) , R(y) = 1 RV . (5. 16) 

/? y - X - BY i 

To obtain the loop equation we start with the identity 

W(x) - R{y) = {y- x)fl(x, y) - BU^x, y), (5.17) 
where we denoted by R(y) the one-point function with Dirichlet boundary condition, 

J?(x) = i(trR( 1 /)) J (5.18) 

and introduced the auxiliary function 

n 1 (x,y) = UtTY 1 W(x)K(y)). (5.19) 

P 

The function Q,\ satisfies the identity 

n 1 (x,y) = B <p — , (5.20) 

which follows from (A.3). After symmetryzing with respect to x we get 

Qx(x, y) + n^-x, y)=B Q(x, y)Q(-x, y). (5.21) 
From here we obtain a quadratic functional equation for the correlator f2: 

(x - y)n(x, y)-(x + y)Sl(-x, y) + W(x) + W(-x) + B 2 Q(x, y)Q(-x, y) = 2R{y). (5.22) 
The linear term can be eliminated by a shift 

G(x,y)=Bn(x,y)-^-. (5.23) 

The function G satisfies 

2 2 

G(x, y)G(-x, y) = —W(x) - W(-x) + 2R(y) - . (5.24) 

This equation is to be compared with the loop equation (5.14) for Ord/DJS boundary conditions, with 

nm = 1 and A (2) = 0: 

J-AW(.x,y)BW(-x,y) = - W(x) - W{-x) + H(y) - ^ . (5.25) 
A(i) A(i) 

These two equations coincide in the limit B, A(i) — > oo. To see this it is sufficient to notice that 
in the limit A(i) — > oo we have the relation = yA^\ The exact relation between A^ (with 
nW = 1, A(2) = 0) and G in this limit follows from the definitions of Dq and Q: 

A^(x,y') = B G ^ v) - G ^- y \ \ W =B\ y' = y 2 (B - oo). (5.26) 
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6 Scaling limit 



In this section we will study the continuum limit of the solution, in which the sum over lattices is 
dominated by those with diverging area and boundary length. The continuum limit is achieved when 
the couplings x, y and p, are tuned close to their critical values. 

Once the bulk coupling constants are set to their critical values, we will look for the critical line 
in the space of the boundary couplings y, A(i) and A(2). After the shift (4.7) the bondary cosmological 
constant x has its critical value at x = 0, while the critical value of y in general depends on the values 
of A(i) and A<2) . 

6.1 Scaling limit of the disk one-point function 

Here we recall the derivation of the continuum limit of the one-point function W(x) from the func- 
tional equation (5.3). Even if the result is well known, we find useful to explain how it is obtained in 
order to set up the logic of our approach to the solution of the functional equations (5.14). 

In the limit x — > 0, the boundary length |<9T| of the planar graphs in (3.7) becomes critical. The 
quadratic functional equation (5.3) becomes singular at x — > when the coefficient A on the r.h.s. 
vanishes. This determines the critical value of the cosmological constant j2, for which the volume |T| 
of a typical planar graph diverges. The condition that the coefficient B of the linear term vanishes 
determines the critical value of the temperature coupling T = T c for which the length of the loops 
diverges: 

r c = l + ^€[l,2]. (6.1) 

Near the critical temperature the coefficient B is proportional to T — T c . 

We rescale x — > ex, where e is a small cutoff parameter with dimension of length, and define the 
renormalized coupling constants as 

A - Mc ~ e 2 /i , T - T c ~ e 26 t (6.2) 

and write (5.1) as 

W(x) = W reg + e 1+e w(x) . (6.3) 

The renormalized bulk and boundary cosmological constants are coupled respectively to the renor- 
malized area A and boundary length i of the graph Y defined as 

A = e 2 |r|, e = e\dr\. (6.4) 

In the following we define the dimensions of the scaling observables by the way they scale with x. We 
will say that the quantity / has dimension d if the ratio f /x d is invariant with respect to rescalings. In 
this case we write [/] = d. The continuous quantities introduced until now have scaling dimensions 

[x] = 1, \pt] = 2, [t] = 29, [to] =1 + 0. (6.5) 

The scaling resolvent w(x) is a function with a cut on the negative axis in the x -plane. It can 
be obtained from the general solution found in [32] by taking the limit in which the cut extends to 
the semi-infinite interval [—00, —M]. To determine M as a function of \i and t one has to solve a 
system of difficult transcendental equations. A simpler indirect method was given in [33]. We begin 
by noticing that the term Cx A in (5.3) drops out because it vanishes faster than the other terms when 



15 



x — ► 0, and B = B\t, where B\ depends only on n. Introducing a hyperbolic map which resolves 
the branch point at x = —M, 

x = Mcoshr , (6.6) 
we obtain a quadratic functional equation for the entire function w(t) = w[x(t)]: 

w 2 (t + +w 2 (t) + nw(r + i-K)w(r) = A + BitM 2 cosh 2 r . (6.7) 

This equation does not depend on the cutoff e, which justifies the definition of the renormalized 
thermal coupling in (6.2). Then the unique solution of this equation is, up a factor which depends on 
the normalization of t, 

W ( T ) = M 1+e cosh(l + 9)t + tM 1 ' 9 cosh(l - 9)t . (6.8) 

One finds B 1 = 4sin 2 (7r60 and A = sm 2 (7r9)(M 1+e - tAl 1 " 9 ) 2 for this solution. 

The function M = M([i, t) can be evaluated using the fact that the derivative d IM W(x) depends 
on /x and t only through M. As a consequence, in the derivative of the solution in /i at fixed x, 

daw = —MdnM ( (1 + 6)M e - (1 - 9)tM~ e ) !E^I (6 .9) 
V / smhr 

the factor in front of the hyperbolic function must be proportional to M : 

B^mUi + 6)M 9 - (1 - 6)tM^ ~ M e ~ l . 
Integrating with respect to fi one finds, for certain normalization of \i, 

H= {I + 9)M 2 - tM 2 ~ 2e . (6.10) 

To summarize, the disk bulk and the boundary one-point functions with ordinary boundary condi- 
tion, —d^U and —d x ll, are given in the continuum limit in the following parametric form: 

-d x U\^ = Af 1+e cosh((l + 6)t) +tM 1 ~ e cosh(l-9)T, (6.11) 
-d^U\ x ~ M e cosh 9t, 
x = M coshr, 

with the function M(/i,t) determined from the transcendental equation (6.10). The expression for 
d^U was obtained by integrating (6.9). 

The function M(t, //) plays an important role in the solution. Its physical meaning can be revealed 
by taking the limit x — ► oo of the bulk one-point function —d^U (x). Since x is coupled to the length 
of the boundary, in the limit of large x the boundary shrinks and the result is the partition function of 
the 0(n) field on a sphere with two punctures, the susceptibility u(fi,t). Expanding at x — > oo we 
find 

- dyU ~ x e - M 2e x~ e + lower powers of x (6. 12) 

(the numerical coefficients are omitted). We conclude that the string susceptibility is given, up to a 
normalization, by 

u = M 26 . (6.13) 
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The normalization of u can be absorbed in the definition of the string coupling constant g s ~ \j [3. 
Thus the transcendental equation (6. 10) for M gives the equation of state of the loop gas on the sphere, 

(1 + 9)vk - tvnr = fx . (6.14) 

The equation of state (6.14) has three singular points at which the three-point function of the 
identity operator d^u diverges. The three points correspond to the critical phases of the loop gas on 
the sphere. At the critical point t = the susceptibility scales as u ~ fi e . This is the dilute phase 
of the loop gas, in which the loops are critical, but occupy an insignificant part of the lattice volume. 
The dense phase is reached when t/x e — > — oo. In the dense phase the loops remain critical but 

1-9 

occupy almost all the lattice and the susceptibility has different scaling, u ~ li « . The scaling of the 
susceptibility in the dilute and in the dense phases match with the values (2.4) and (2.6) of the central 
charge of the corresponding matter CFTs. Considered on the interval — oo < t < 0, the equation of 
state (6. 14) describes the massless thermal flow [34] relating the dilute and the dense phases. 
At the third critical point d^M becomes singular but M itself remains finite. It is given by 

t e = ^ Mf > 0, li c = M 2 C < 0. (6.15) 

Around this critical point /j — Li c ~ (M — M c ) 2 + • • • , hence the scaling of the susceptibility is that 
of pure gravity, u ~ (/i — /U c ) 1//2 . 



6.2 The phase diagram for the DJS boundary condition 

We found the scaling limit of the one-point function (4.11) as a function of the renormalized bulk 
couplings, li and t, and the coupling x characterizing the ordinary boundary. Now, analyzing the 
loop equation (5.14) for the two-point functions, we will look for the possible scaling limits for the 
couplings y, A(i) and A (2) , characterizing the DJS boundary. 

As in the previous subsection, we will write down the conditions that the regular parts of the 
source terms vanish. Let us introduce the isotropic coupling A and the anisotropic coupling A as 

A ( i) = A + \ A, A( 2 ) = A - |A (6.16) 

and substitute (5.1) in the r.h.s. of (5.14). We obtain 

= Co + c\x + C2X 2 — A(w + nWw) , 

= c + cix + c 2 x 2 + A(w + nWw) , (6.17) 
where the coefficients Co and c\ are functions of A and A: 

n-i (77(1) — 77 ( 2 ) ) 

co = (A - \^)( 92 H + 93^1) + 2Ag + y - Aift - A - - ' (6.18) 

Cl = (A'-iA'ta^-l-^ + A ^J^ . (6.19) 

/ n (l) _ n (2) \ 

For generic values of the couplings y, A and A, the coefficient co is non-vanishing. The condition 
cq = determines the critical value y c where the length of the DJS boundary diverges. 6 Once the 



6 Indeed, the term co is the dominant term when x — > 0. For Co 7^ the solution for and B (a ' 1 in (5.14) is given by 
linear functions of x and w. Such a solution describes the situation when the length of the DJS boundary is small and the 
two-point function degenerates to a one-point function. 
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boundary cosmological constant is tuned to its critical value, the condition c\ = determines the 
critical lines in the space of the couplings A(i) and A<2), where the DJS boundary condition becomes 
conformal. The two equations 

c (y,A,A) = 0, Cl (y,A,A)=0 (ji = t = 0) (6.21) 

define a one-dimensional critical submanifold in the space of the boundary couplings {y, m, A}: 

A = A*(A). (6.22) 

Obviously and A^ cannot be simultaneously zero. Therefore the curve (6.22) consists of two 
branches, which correspond to different conformal DJS boundary conditions, the lines of anisotropic 
special transitions AS^ and AS( 2 )'- 

AS {1) : A {1) = 0, A {2) + 0; 

AS {2) : A (2) = 0, A {1) + 0. (6.23) 

The branch AS^ corresponds to A > 0, while the branch AS (2) corresponds to A < 0. Consider 
the behavior of correlators on the two branches of the critical line. By the definition (5.15), the 
two correlators and are related by 

n( 2 ) n W 
DP = -J*, D f) = _5b . (6.24) 
1 - AD {2 > 1 + AD£> 

On the branch ASm the correlator diverges while remines finite, and vice versa. Assume 
that A(i) and A (2 ) are positive. Then Dq*' are both positive by construction. If A = A(i) — A(2) > 0, 



then the coefficients of the geometric series 

OO , 

D^ = Y^^ k (D^) (6.25) 



k=0 

are all positive and diverges, while — ► 1/A. Thus the branch AS^, where diverges, 
is associated with A > 0. On this branch the probability that the loops of color (1) touch the DJS 
boundary is critically enhanced. In the correlator Djp, the ordinary boundary behaves as a loop of 
type (1) and can touch the DJS boundary. The geometrical progression (6.25) reflects the possibility 
of any number of such events, each contributing a factor A. Conversely, the ordinary boundary for the 

(2) 

correlator D behaves as a loop of type (2), since such loops almost never touch the DJS boundary. 
On the branch AS( 2 ) the situation is reversed. 

It is not possible to solve explicitly the conditions of criticality (6.21) without extra information, 
because they contain two unknown functions of the three couplings, H and H±. Nevertheless, the 
qualitative picture can be reconstructed. 

First let us notice that the form of the critical curve can be evaluated in the particular cases nW = n 
and nW = 0. In the first case n( 2 > = and the correlation functions do not depend on A(2) and so the 
coefficients c\ and c 2 depend on A and A through the combination A(i) = A + A/2. Similarly one 
considers the case nW = 0. The phase diagram in these two cases represents an infinite straight line 
separating the ordinary and the extraordinary transitions: 

f A c - A/2 if nW = n, 
A (A) = < 1 (6.26) 

IA C + A/2 if nW = 0. 
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X 



(2) 



Figure 5: Phase diagram in the rotated (A(i>, A<2>) plane for n < 1 and n< x > > n (2) . The ordinary and the 
extraordinary phases are separated by a line of anisotropic special transitions, which consists of two branches, 
ASn) and AS^)- The two extraordinary phases, Ext^ and Ext^), are separated by the isotropic line A = 

The critical line crosses the axis A = at the special point A = A c . The value of A c can be evaluated 
by solving (6.21) for n* 1 ) = n and A(2) = 0. The result is 

K = v ' (6 + " )(2 -" ) . (6.27) 

1 — n 

For general nW g [0, n] we can determine three points of the critical curve: 

1 — TX 1 — Tl 

{X W ,\ {2) } = { T ——\ C ,0}, {0, T -^A C }, {A C ,A C }. (6.28) 

In the two limiting cases considered above the critical line, given by equation (6.26), crosses the 
anisotropic line without forming a cusp. Is this the case in general? Let us consider the vicinity of the 
special point (A, A) = (A c , 0). In the vicinity of the special point a new scaling behavior occurs. In 
this regime the term x 2 in (6.17) cannot be neglected. The requirement that all terms in (6.17) have 
the same dimension determines the scaling of the anisotropic coupling A: 

[A] = 1 - 6. (6.29) 

In order to determine the form of the critical curve near the special point, we return to the equations 
(6.21) and consider the behavior near the special point of the unknown functions H(y) = 4 ( trH(y) ) 
and H\(y) = ^ ( trXH(y) ), which depend implicitly on A and A. For t = \x = 0, these functions can 
be decomposed, just as the one-point function with ordinary boundary condition, W(x), into regular 
and a singular parts: 

H(y)=H^(y) + h(y), H^y) = H?*(y) + est ■ h(y) . (6.30) 
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On the critical curve A = A* (A) the singular part of H vanishes and the coefficient c\ given by eq. 
(6.19) can be Taylor expanded in A — A c and A: 

ci(m, A) = A x (A - A c ) + B 1 (n« - n( 2 >)A + A 2 (A - A c ) 2 + B 2 A 2 + • • • = 0. (6.31) 

If A\ ^ 0, the critical curve is given by a regular function of A and A, the critical curve is a continuous 
line which crosses the real axis at A = A c without forming a cusp. This form of the curve differs from 
the predictions of [2] and [3], where a cusp-like form is predicted by a scaling argument. We will see 
later that the fact that the critical curve is analytic at the special point does not contradict the scaling 
(6.29). 



6.3 Scaling limit of the functional equation for the disk two-point function 
6.3.1 The scaling limit for A^O 

Consider first the case when the anisotropic coupling A is finite and assume that we are on the branch 
AS(i^ where A > 0. Then the x 2 term on the r.h.s. can be neglected, because it is subdominant 
compared to w ~ x 1+d . The scaling limit corresponds to the vicinity of the critical submanifold 
where the two coefficients cq and c\ scale respectively as x 1+e and x e . 

We are now going to find the scaling limit of the loop equations (5.14). In the scaling limit we can 
retain only the singular parts of the correlators which we denote by 4^'^ (L = 0, 1, ...). 

We define the functions d^ and d^f 1 by 



,(i) _ w« n (2) L_ ^ I _l 1_ J2) 

'o 



u Q — a Q , u Q — ~ -t- ^ 2 a Q . {o.jz) 



Then the relation (6.24) implies 



4 1} 4 2) = -1. (6.33) 



We define in general 4^'^ as the singular part of , with the normalization chosen so that the 

recurrence equation (5.7) holds for any L > 0: 

4+i = w-kdf, 

4li = ™*4 2) - ( 6 - 34 ) 

On the branch ASru the first of the two equations (5.14) becomes singular, since A^ vanishes while 
A^ remains finite. We write this equation in terms of d^ and d^ using that 

AU = \ m d£\ B W = (6 .35) 

A(i) 

We get 

4 2) 4 } +w + nWw = hb -t B x , (A > 0) (6.36) 

where \i b and t b are defined by 

-r = l*B, A = ~ 



20 



Once the solution of (6.37) is known, all two-point functions rife/ can be computed by using the 
recurrence equations (6.34). 



The scaling limit near the branch ASf 2 ) (A < 0) is obtained by using the symmetry raf 1 ) <-> 
n< 2 > , A <-> —A. In this case one obtains another equation 

d { Q ] df ] +w + nV)w = fj, B -t B x . (A<0) (6.38) 

Note that the relation (6.33) is true on both branches of the critical line. 

The map {y, A} — > {^b^b} defined by (6.18), (6.19) and (6.37) represents a coordinate change 
in the space of couplings which diagonalizes the scaling transformation. The coupling \i B is the 
renormalized boundary cosmological constant for the DJS boundary. 7 The coupling ig is the renor- 
malized boundary matter coupling, which defines the DJS boundary condition. The dimensions of 
these couplings are 

\Hb] = 1 + 0, [t B ] = 9. (6.39) 

Once we choose y so that [lb = 0, the condition t B = gives the critical curve where the anisotropic 
special transitions take place. If the function t B is regular near the critical line A = A* (A), then it can 
be replaced by the linear approximation 

t B ~ A* (A) - A. (6.40) 

The deformations in the directions t B and A, are driven by some Liouville dressed boundary 
operators Of and 0\. Knowing the dimensions of t B and A, we can determine the Kac labels of 
these operators with the help of the KPZ formula. The general rule for evaluating the Kac labels in 
2D gravity with matter central charge (2.4) is the following. If a coupling constant has dimension a, 
then the corresponding operator has Kac labels (r, s) determined by 

/ e (i + e)r-s\ 

a = « riS = min l + -±^ ± . (6.41) 



The details of the identification are given in Appendix B. We find 

tn 

Near the special point we have 



Ol = 0* Ol = Of*. (6.42) 



A~t^, 0=-iL = ^<l. (6.43) 

Since t B and A scale differently, there is no contradiction between the scaling (6.43) and the analyticity 
of the critical curve near the special point. 



7 More precisely, it is a combination of the boundary coupling constant and the disk one-point function with DJS bound- 
ary. What is important for us is that the condition /is = fixes the critical value of the bare DJS cosmological constant y. 
At fiB = 0, the length of the DJS boundary diverges. 
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6.3.2 The scaling limit in the isotropic direction (A = 0) 

Along the isotropic line A = the two functional equations (5.14) degenerate into a single equation 
for the correlator Dq: 

As , D „-^- y-*+ff^-^ . (,44, 

y — x + A [W + H - V) 

In order to evaluate D\ = = Df \ we can consider the linear order in A. It is however easier to 
use the fact that Dq and D\ do not depend on the splitting n = nW + n< 2 ). Furthermore, if we choose 
nW = n and A(i) = A, the observables do not depend on A(2) , which can be chosen to be zero. Taking 
= n, A(i) = m and = A(2) = 0, we obtain from (5.14) 

y - x + \{H - W - nW) 



# (1) U=n = A ( XD ^ +nW)+x-y= y - ^D - 1 " (6 " 45) 

From these expressions it is clear how the scaling of the singular parts of Dq and D\ , which we 
denote respectively by do and d\, change when we go from A = to A = A c . When A = we 
have H(y) = W(y) and Dq is the disk partition function with ordinary boundary conditions and two 
marked points on it, eq. (4.17). When A = A c and y = y c , 

2 
7" 

w g 3 x< 



w(w + nw) O o 

d = g 3 —~x 1 - e , d x = \ 2 > ~x 29 (X = \ c ,y = y c ). (6.46) 



6.3.3 Dirichlet versus DJS 

Now we will focus on the special case nt 1 ) = 1 and compare the scaling behavior with that for the 
Dirichlet boundary conditions. The critical behavior of the two-point correlator in both cases is the 
same, but the boundary coupling constants correspond to different boundary operators. 

Consider the functional equation (5.25) for the correlator with Ord/DJS boundary conditions when 
nW = 1. The critical value of A(i) is infinite in this case, see equation (6.28). The scaling limit of 
(5.25) is 

% {x,y)d\ l (-x,y) = - w(x) - w(-x) + \i B - t — • (6.47) 

A(i) 

The last term remains finite if A(i) tends to infinity as x~ 6 . The scaling boundary coupling can be 
identified as tg = 1/A(i) and equation (6.47) takes the general form (6.36). What is remarkable here 
is that the boundary temperature constant need not to be tuned. Equation (6.47) holds for any value 
of A(i). On the other hand, when ts = 1/A(i) is small, the last term describes the perturbation of 
the ASfu boundary condition by the thermal operator 0\ .3 with ai 3 = 9. When A(i) is small, the 
last term accounts for the perturbation of the ordinary boundary condition by the two-leg boundary 
operator O31 with 031 = —6, whose matter component is an is irrelevant operator. 

Now let us take the scaling limit of the quadratic functional equation (5.24) for the correlator with 
Ord/Dir boundary conditions. At x = the equation (5.24) becomes algebraic. The critical value 
y = y c , where the solution develops a square root singularity, is determined by 

2W(0) - 2R(y c ) + y 2 c /B 2 = 0. 

We can write equation (5.24) as 

G(x, y)G(-x, y) = fi B - w (x) - w(-x) - , (6.48) 
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where 

2 2 

fi B = 2w(0)+2[R(y)-R(y c )] + ^^. (6.49) 

For any finite value of B, the scaling limit of this equation is 

G(x, y)G(-x, y) = fi B - w{x) - w(-x) . (6.50) 

The x 2 term survives only if B vanishes as x~ 

[B] = (1 - 9)/2 = a 2 ,i. (6.51) 

This is the expected answer, because \ is the one-leg boundary operator which creates an open line 
stalling at the boundary. We conclude that the Dirichlet and the DJS boundary conditions have the 
same scaling limit, but in the first case the boundary coupling A corresponds to a relevant perturba- 
tion and it is sufficient give it any finite value, while in the second case the boundary coupling A(i) 
corresponds to an irrelevant perturbation and therefore must be infinitely strong. 

7 Spectrum of the boundary operators 

Let us denote by and the scaling dimensions respectively of a£ and : 

4 Q) = [4 a) ] (L>0, a = 1,2). (7.1) 
The recurrence equations (6.34) tell us that the dimensions grow linearly with L: 

4 Q) = L[w] + <4 a) . (7.2) 

These relations make sense in the dilute phase, where [w] = 1 + 9, as well as in the dense phase, 
where [w] = 1 — 9. In addition, by (6.33) we have 

a { Q ] + af> = 0. (7.3) 

Thus all scaling dimensions are expressed in terms of off . 

Let us evaluate ao for the branch AS^ of the critical line. We thus assume that A is finite and 
positive sufficiently far from the isotropic special point. Take /i = /xb = ts = and write the shift 
equations which follow from (6.36), 

d^\e lw x) _ w(e~ t7T x) + nWw(x) 
(1) ' 4\e-™ X ) ~ w(e™x)+nm W (x) ' 

The one-point function (6.8) behaves as w ~ x l+e in the dilute phase (t = 0) and asw~ x 1 " 61 in the 
dense phase (t — > — oo). In both cases the r.h.s. is just aphase factor. Since all the couplings except for 

x have been turned off, (x) should be a simple power function of x. Substituting (x) ~ x a L 
in (7.4), we find 

sinvrfan 1 ^ ± 9) 

n w = (1) , (+ for d i lute > - for dense ) • ( 7 -5) 

sin -rag 
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This equation determines the exponent ckq up to an integer. In the parametrization (2.15) we have 

= —Or + jdii in the dilute phase and = Or + jden in the dense phase. 
The integers and j^ en can be fixed by additional restrictions on the exponents. Let us assume 
that A(i) and A(2) are non-negative, n > and the boundary parameter r is in the 'physical' interval 
1 < r < 1/9—1, where both nW and n< 2 ) are non-negative. These assumptions guarantee that 
the Boltzmann weights are positive and the loop expansion of the observables has good statistical 
meaning. Since all loop configurations that enter in the loop expansion of the one-point function 
W{x) are present in the loop expansions of A {a) andSH the singularity of these observables when 
A — > A* ( A) must not be weaker than that of W. In other words, the scaling dimensions of ~ A^ 



and dp ~ must not be larger than the scaling dimension of the one-point function 



w: 



a [ Q ] < [w], aS 1} < [w]. (7.6) 

Since ckq + a£ = [w], this also means that c^f 1 and are non-negative. Taking into account that 
[w] = 1 ± 6 in the dilute/dense phase, we get the bound 

- (1 ± 9) < a$ < (+ for dilute, - for dense) . (7.7) 

This bound determines jdii = and j'den = —1- As a consequence, on the branch AS^ of the critical 
line the dimensions afip = [d^] in the dilute and in the dense phases are given by 

AS . = L(l + 9) - 9r, af = L(l + 0) + Or (dilute phase) 

(1) ' a { l ] = L(l - 9) +9r - 1, = L(l - 0) - Or + 1 (dense phase). 

Note that the results for the dense phase are valid not only in the vicinity of the critical line AS/]\, but 
in the whole half -plane A > 0. 

By the symmetry (1) <-> (2), the exponents a.^ on the branch ASn-\ and the exponents on 
the branch AS^) should be related by nW <-> n< 2 ), or equivalently r <-> 1/0 



r: 



AS . = L(l + 0) - Or + 1, af = L(l + 0) + Or - 1, (dilute phase) 

(2) ' o4 1} = L(l -0) + Or, a { l ] = L(l - 0) - Or (dense phase). 

The scaling exponents of the two-point functions of the 0(n) model coupled to 2D gravity al- 
low, through the KPZ formula [23, 24], to determine the conformal weights of the matter boundary 
operators. In the dilute phase, where the Kac parametrization is given by (2.5), the correspondence 
between the scaling dimension a of a boundary two-point correlator and the conformal weight h rjS of 
the corresponding matter boundary field is given by 

a = (1 + 0)r - s -> h = h TjS (dilute phase). (7.10) 

In the dense phase, where the Kac labels are defined by (2.7), one obtains, taking into account that the 
identity boundary operator for the ordinary boundary condition has 'wrong' dressing, 

a = r — s(l — 0) — > h = h T:S (dense phase). (7.11) 

From (7.10) and (7.11) we determine the scaling dimensions of the L-leg boundary operators (3.9): 

AS . sP^O?_ L>r , 4 2) ^O r B +L , r (dilute phase) 

4 1} ^ ®r-i,r-L, 4 2) - Of-i,r+L (dense phase). 
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AS . - Of_ w S<jf> ^ 0*. Ltr+1 (dilute phase) 

(2) ' 4 1} - ^r-L- 4 2) - ^r+L (dense phase). 

These conformal weights are in accord with the results of [7], [13], [14], [3]. We remind that the 
scaling dimensions are determined up to a symmetry of the Kac parametrization: 

hr t s = h- r - s , hr,s-i = K+i/e,s+i/e (dilute phase) (7.14) 
h r , s = h- r - s , h r+ i )S = h r+1 /Q tS+1 /0 (dense phase). (7.15) 

Comparing the scaling dimensions in the dilute and in the dense phase, we see that the bulk thermal 
flow t O13 transforms the boundary operator O r s in the dilute phase into the boundary operator 
O s -i >r in the dense phase. For the rational points 9 = 1/p, our results for the endpoints of the bulk 
thermal flow driven by the operator 1 0\ $ match the perturbative calculations performed recently in 
[15]. 

In the O(n) model the boundary parameter r is continuous and we can explore the limit r — > 1, 
in which the BCC operator O r>r carries the same Kac labels as the identity operator. Let us call 
this operator Of 1 . The bulk thermal flow transforms the operators Of 1 and Of 1 into two different 
boundary operators in the dense phase. Hence there are at least two distinct boundary operators with 
Kac labels (1,1): the identity operator and the limit r — > 1 of the operator O l r 



r,r- 



8 Solution of the loop equations in the scaling limit 

We are going to study two particular cases where the analytic solution of the functional equations 
(6.36) and (6.38) is accessible. First we will evaluate the two-point function on the two branches 
ASt]\ and ASr 2 \ of the critical line, where the 0(n) field is conformal invariant both in the bulk and 
on the boundary. In this case t = t B = and the boundary two-point function is that of Liouville 
gravity. The three couplings are introduced by the world sheet action of Liouville gravity with matter 
central charge (2.4), which we write symbolically as 

5Liouv = cSfree + / ^1,1+ J X Of^ + j fi B 0? >v (8.1) 
bulk Ord. boundary DJS boundary 

Since the perturbing operators in this case are Liouville primary fields, 



~ e 2 ^, Of 



the two-point function is given by the product of matter and Liouville two-point functions. Up to a 
numerical factor, the solution as a function of fi and \l# must be given by the boundary two-point 
function in Liouville theory [35]. We will see that indeed the functional equation (6.36) is identical 
to a functional equation obtained in [35] using the operator product expansion in boundary Liouville 
theory. 

In the second case we are able to solve, we take fi = = and non-zero matter couplings t and 
ts- This case is more interesting, because it is not described by the standard Liouville gravity. The 
corresponding world sheet action is symbolically written as 

5 = 5 fre e+ J tO h3 + J xOf A + j t B Of 3 . (8.2) 

bulk Ord. boundary DJS boundary 



25 



The worldsheet theory described by this action is more complicated than Liouville gravity, because it 
does not enjoy the factorization properties of the latter. The boundary two-point correlator does not 
factorize, for finite t and/or t B , into a product of matter and Liouville correlators, as is the case for the 
action (8.1). This is because the perturbing operators and Of% have both matter and Liouville 
components: 

01,3 ~ *i,3 e 2b ^t 0* ~ $g 3 e* 1 -*)*. 

Let us mention that the theory of random surfaces described by the action (8.2) has no obvious direct 
microscopic realization. Our solution interpolates between the two-point functions for the dilute (t = 
0) and the dense (t — ► — oo) phases of the loop gas, on one hand, and between the anisotropic special 
(t B = 0) and ordinary /extraordinary boundary conditions (t B — > ±oo), on the other hand. 

8.1 Solution for t = 0, t B = 

In the dilute phase (t = 0) the solution (6.6) - (6.8) for the boundary one-point function takes the form 
x = M cosh r, w(x) = M 1+e cosh(l + 9)t. (8.3) 
Then the loop equations (6.36) become a shift equation 

d o ] ( T ) d i\ T ± i7T ) + ^ocosh [(1 + 0)t ± »7r(l - r)6] = fj, B - t B M coshr, (8.4) 
where we introduced the constant 

Wo = M ^^0_ (85) 
sin irrv 

At the point t B = 0, where the DJS boundary condition is conformal, this equation can be solved 
explicitly. After a shift r — ► r =F in we write it, using (6.33), as 

4 1} (r) = [too cosh [(1 + 0)r ± i7rr0] -/i B ] d { Q l) (t ± iir) . (8.6) 

If we parametrize fi B in terms of a new variable a as 

//B = wo cosh(l + 0)a, (8.7) 

the loop equation turns out to be identical to the functional identity for the boundary Liouville two- 
point function [35], which we recall in Appendix B. In the Liouville gravity framework, r and a 
parametrize the FZZT branes corresponding to the ordinary and anisotropic special boundary condi- 
tions. 

The loop equations for the dense phase (t — > — oo), are given by (8.6) with 6 sign-flipped. This 
equation describes the only scaling limit in the dense phase. The term with t B is absent in the dense 
phase, because it has dimension 1 + 9, while the other terms have dimension 1 — 9. In this case, the 
loop equation gets identical to the functional identity for the Liouville boundary two-point function if 
we parametrize fi B as 

Mb = cosh (l — 9) a. (8.8) 
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8.2 Solution for = 0, /x B = 

Here we solve the loop equation (6.36) in the scaling limit with fi = \ig = but keeping x, t and ts 
finite. Let us first find the expression for the one-point function w(x) for \i = 0. The equation (6.10) 
has in this case two solutions, M = and M = (1 + 9)~ ia?. One can see [33] that the first solution 
is valid for £ < 0, while the second one is valid for t > 0. Therefore when // = and t < 0, the 
solution (6.6)-(6.8) takes the following simple form: 

= x 1 ^ + tx 1 ' 9 (t < 0). (8.9) 

Introduce the following exponential parametrization of x,t,tB in terms of r, 7, 7: 

s = e r , i = - e 2 ^, t B = -2w e^smh(j9), w = ^t^t - (8.10) 

sm(7rr#) 

In terms of the new variables, equation (6.36) with (jlb = A 4 = acquires the form 

diVyd^Cr ± iTr) = «;„ (-2e* sinh^K + e (i+^±«^ _ e 2*7 e (i+tf)r±i«r^ 

= 4w f+fi cosh g(T " 7 + 1* sinh 9{T ~ 7 " 1 ± mr) . (8.11) 



Taking the logarithm of both sides we obtain a linear difference equation, which can be solved explic- 
itly. The solution is given by 

4\r) = i e -i-^-^ + iF_ r (T-7 + 7)^- r (r-7-7) ; 
d?(r) = -el^(H^Hiy 1 _ r (r-7 + 7)V 1+| _ r (r-7-7), 

where the function y r (r) is defined by 



, „ , , def 1 f du 
lo g y r (r) = - - 



sinh(vrrLij) rl 



sinh(7ro;) sinh ttlo 



(8.13) 



The properties of the function V r (r) are listed in Appendix C. 

The solution (8.12) reproduces correctly the scaling exponents (7.8) and it is unique, assuming 

that the correlators and d^ have no poles as functions of x. Near the branch AS(2), the functions 

(2) (2) 

dy and d\ are given by the same expressions (8.12), but with r replaced by 1/9 — r. 
8.3 Analysis of the solution 

To explore the scaling regimes of the solution (8. 12) we use the expansion (C.2) and return to the 
original variables, 



f = x, = (-*)*, e^=iT^ + ff -tj . (8.14) 

Let us define the function V(x) by V r {r) = V r {e T ). The large x expansion of V r goes, according to 
(C.5), as 

* . . ^o/o ( sin7rr 1 sin7r7"# fl \ 

V r (x) = x re ' 2 1 + x- 1 + X~ + ...). (8.15) 
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The expansion at small x follows from the symmetry V(x) = V(l/x). Written in terms of the original 
variables, the scaling solution near the branch ASm\ takes the form 



4\ X ) = — 



w 



t B /2 + Jt 2 B /4 - t 



x V- 



x { t B /2 + Jt 2 /4-t 



4 1} (x) 



l-r 
Xt 2 



-t B /2 + Jt 2 B /4 - t 



x ( t B /2 + Jt 2 B /4 - t 



Vi 

H 

21) 



-t B /2 + Jt 2 B /4 - t 



(8.16) 



-t B /2 + Jt 2 B /4 - t 



The critical regimes of this solution are associated with the limits t 
of the bulk and the boundary temperature couplings. 



.(8.17) 

0, — oo and t B — *■ 0, ±oo 



( i) Dilute phase, anisotropic special transitions 

This critical regime is achieved when both t and t B are small. Using the asymptotics (8.15), we 
find that in the limit (t B — > 0, t — > —0) the expressions (8.16)-(8.17) reproduce the correct scaling 
exponents (7.8) in the dilute phase: 



AS, 



(i) 



d« ~ x 



d W ^ x i+e-er 



B 



0, t 



-0). 



(8.18) 



The regime AS( 2 ) i s obtained by replacing (1) — ► (2), r 



r. 



( ii) Dilute phase, ordinary transition 

At t — ► — 0, the leading behavior of = l/d ^ and d± for large t B is (we omitted all numerical 



coefficients) 

J 2 ) 



7(1) 



t B r (l + t B ~iix + t B l X 6 + ... 

tB 1 ^ (x + t B -^X 2 + i^ 1 + t^ 2 X 20+1 + 

7(2) 



(t -» -0, t B -» +oo). (8.19) 



In the expansion for c?q , the first singular term, x , is the singular part of Do with ordinary/ordinary 
boundary conditions. In the expansion for d± , the first singular term, x l+e , is the one-point function 
w, while the next term, x 2e+1 , is the singular part of the boundary two-point function Dj^ with 
ordinary /ordinary boundary conditions. 



( Hi) Dilute phase, extraordinary transition 

Now we write the asymptotics of (8. 16)-(8. 17) in the opposite limit, t 



d, 



(2) 


7(1) 



t B 



X + t B ~ex £ + t B X a+1 + t R X 



1 ryO + 1 



-2 ^29+1 



+ ... 



d\ X) ~ t B * +1 - r (l + t B -iix + t B 1 x e + ... 



> —0, and t B 
(t = 0,t B - 



-oo: 



-oo 



This asymptotics reflects the symmetry of the solution (8.12) which maps 



1 + 1/0- r, 



7(2) 



7(D 



(8.20) 
(8.21) 
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which is also a symmetry of the loop equations (5.14). In the limit of large and negative ts, the 
function behaves as the singular part of the correlator D^' with ordinary/ordinary boundary con- 
ditions, while behaves as the singular part of the correlator Dq with ordinary/ordinary boundary 
conditions. 

The asymptotics of the solution at ts — ► ±00 confirms the qualitative picture proposed in [3] and 
explained in the Introduction. When ts is large and positive, the loops avoid the boundary and we 
have the ordinary boundary condition. In the opposite limit, ts — ► —00, the DJS boundary tends to be 
coated by loop(s). Therefore the typical loop configurations for D± in the limit ts — ► — 00 will look 
like those of in the ordinary phase, because the open line connecting the two boundary-changing 
points will be adsorbed by the DJS boundary. Conversely, the typical loop configurations for Dq 
will look like those of D± in the ordinary phase, because free part of the loop that wraps the DJS 
boundary will behave as an open line connecting the two boundary-changing points. 

We saw that the solution reproduces the qualitative phase diagram for the dilute phase, shown in 
Fig. 5. Now let us try to reconstruct the phase diagram in the dense phase. 

(iv) Dense phase, anisotropic special transitions 

For any finite value of ts, the dense phase is obtained in the limit t — > — 00. The asymptotics of 

(8.16) -(8.17) in this limit does not depend on ts'- 

4 1} ~ x r6 - x {-t)^~ r , df 5 ~ aT (1 ~ r)9 {-tf- r+ ^ (*->-oo). (8.22) 

This means that in the dense phase the DJS boundary condition is automatically conformal for any 
value of ts- The boundary critical behavior does not change with the isotropic boundary coupling 
ts, but it can depend on the anisotropic coupling A. The solution (8. 16)-(8. 17) holds for any positive 
value of A. For negative A we have another solution, which is obtained by replacing (1) — ► (2) and 
r — > 1/9 — r. Thus in the dense phase there are two possible critical regimes for the DJS boundary, 
one for positive A and the other for negative A, which are analogous to the two anisotropic special 
transitions in the dilute phase. The domains of the two regimes are separated by the isotropic line 
A = 0. 

The above is true when ts is finite. If ts tends to ±00, we can obtain critical regimes with the 
properties of the ordinary and the extraordinary transitions. 

(v) Dense phase, ordinary and extraordinary transitions 

If we expand the solution (8. 16)-(8. 17) for —t S> x 20 and ts ^> —tx~ e , the singular parts of the 
two correlators will be the same as the correlators with ordinary boundary condition on both sides. 
For example, instead of the term t]^x 9 in the expression for d G , we will obtain t~ l tsx~ e . This is 
the singular part of the correlator Dq with ordinary boundary conditions in the dense phase. Further, 
the asymptotics of the solution in the limit —t ^> x 2e and ts <C tx~ e is determined by the symmetry 
(8.21). This critical regime has the properties of the extraordinary transition, in complete analogy with 
the dilute case. We conclude that the ordinary and the extraordinary transitions exist also in the dense 
phase, but they are pushed to ts — ► ±00. 

Finally, let us comment on the possible origin of the square-root singularity of the solution (8.16)- 

(8.17) at t 2 B = At. This singularity appears in the disordered phase, t > 0, which is outside the domain 
of validity of the solution. Nevertheless, one can speculate that this singularity is related to the surface 
transition, which separates the phases with ordered and disordered spins near the DJS boundary. The 
singularity in our solution has two branches, ts = ±2i/t, while in the true solution for t > the 
negative branch should disappear. 
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9 Conclusions 



In this paper we studied the dilute boundary 0(n) model with a class of anisotropic boundary condi- 
tions, using the methods of 2D quantum gravity. The loop gas formulation of the anisotropic boundary 
conditions, proposed by Dubail, Jacobsen and Saleur (DJS ), involves two kinds of loops having fu- 
gacities n' 1 ) and n< 2 ) = n — nW . Besides the bulk temperature, which controls the length of the loops, 
the model involves two boundary coupling constants, which define the interaction of the two kinds of 
loops with the boundary. 

The regime where the DJS boundary condition becomes conformal invariant is named in [3] 
anisotropic special transition. The enhanced symmetry of the model after coupling to gravity sys- 
tem allowed us to solve the model analytically away from the anisotropic special transition. We used 
the solution to explore the deformations away from criticality which are generated by the bulk and 
boundary thermal operators. 

Our main results can be summarized as follows. 

1) We found the phase diagram for the boundary transitions in the dilute phase of the 0(n) model 
with anisotropic boundary interaction. The phase diagram is qualitatively the same as the one obtained 
in [3] and sketched in the Introduction. The critical line consists of two branches placed above and 
below the isotropic line. Near the special point the critical curve is given by the same equation on both 
sides of the isotropic line, which means that the two branches of the critical line meet at the special 
point without forming a cusp. We also demonstrated that the analytic shape of the critical curve does 
not contradict the scaling of the two boundary coupling constants. This contradicts the picture drawn 
in [2] on the basis of scaling arguments, which seems to be supported by the numerical analysis of 
[3]. Of course we do not exclude the possibility that the origin of the discrepancy is in the fluctuations 
of the metric. 

2) From the singular behavior of the boundary two-point functions we obtained the spectrum 
of conformal dimensions of the L-leg boundary operators between ordinary and anisotropic special 
boundary conditions, which is in agreement with [3]. In order to establish the critical exponents we 
used substantially the assumption that nW and n< 2 ) are both non-negative. 

3) We showed that the two-point functions of these operators coincide with the two-point func- 
tions in boundary Liouville theory [35]. The functional equation for the boundary two-point function 
obtained from the Ward identities in the matrix model is identical to the functional equation derived 
by using the OPE in boundary Liouville theory. 

4) The result which we find the most interesting is the expression for the two-point functions 
away from the critical lines. For any finite value of the anisotropic coupling A, the deviation from 
the critical line is measured by the renormalized bulk and the boundary thermal couplings, t and ts- 
Our result, given by eqs. (8. 16)-(8. 17), gives the boundary two-point function in a theory which is 
similar to boundary Liouville gravity, except that the bulk and the boundary Liouville interactions 
are replaced by the Liouville dressed bulk and boundary thermal operators, 1 0^3 and ts Of 3- The 
boundary flow, generated by the boundary operator Of 3 , relates the anisotropic special transition 
with the ordinary and the extraordinary ones. The bulk thermal flow, generated by the operator 0\^, 
relates the dilute and the dense phases of the 0(n) model coupled to gravity. At the critical value of 
the boundary coupling, the bulk flow induces a boundary flow between one DJS boundary condition 
in the dilute phase and another DJS boundary condition in the dense phase. For the rational values of 
the central charge, the boundary conditions associated with the endpoints of the bulk flow match with 
those predicted by the recent study using perturbative RG techniques [15]. 

Here we considered only the boundary two-point functions with ordinary/DJS boundary condi- 
tions. It is not difficult to write the loop equations for the boundary (n + l)-point functions with one 



30 



ordinary and n DJS boundaries. The loop equations for n > 1 will depend not only on the parameters 
characterizing each segment of the boundary, but also on a hierarchy of overlap parameters that define 
the fugacities of loops that touch several boundary segments. The loop equations for the case n = 2 
were studied for the dense phase in [27]. In this case there is one extra parameter, associated with 
the loops that touch both DJS boundaries, which determines the spectrum of the boundary operators 
compatible with the two DJS boundary conditions. In the conformal limit, the loop equations for the 
(1 + n)-point functions should turn to boundary ground ring identities, which, compared to those 
derived in [36] for gaussian matter field, will contain a number of extra contact terms with coefficients 
determined by the overlap parameters. It would be interesting to generalize the calculation of [27] to 
the dilute case and compare with the existing results [37, 38] for the 3-point functions in Liouville 
gravity with non-trivial matter field. 

The method developed in this paper can be generalized in several directions. Our results were 
obtained for the 0(n) model, but they can be easily extended to other loop models, as the dilute ADE 
hight models. It is also clear that the method works for more general cases of anisotropic boundary 
conditions, with the 0(n) invariance broken to O(ni) x • • • x O(nfc). 

Finally, let us mention that there is an open problem in our approach. The loop equations does 
not allow to evaluate the one-point function with DJS boundary conditions, H(y), except in some 
particular cases. There are two possible scaling limits for this function, which correspond to the two 
Liouville dressings of the identity operator with DJS boundary condition, and it can happen that both 
dressings are realized depending on the boundary parameters. This ambiguity does not affect the 
results reported in this paper. 
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A Derivation of the loop equations 

Here we give the derivation of the loop equations which are extensively discussed in this work. 

We first summarize our technique to derive loop equations. The translation invariance of the 
matrix measure implies, for any matrix F made out of X and Y , the identities 

^(a x F) + (tr[(-y / (X)+Y 2 )F]) = 0, (A.l) 

i<0y o F) + <tr[Y o (XF + FX)]) = 0, (A.2) 

where the derivatives with respect to matrices are defined by <9xF = dFij/dXij summed over the 
indices i and j, and are generally given by sums of double traces. Written for the observable G = 
— (XF + FX), the second equation states 8 

<tr(Y a G)> = - £_(a Yo [W(x)GW(-x)]>. (A.3) 

8 Here it is assumed that the eigenvalues of X are all in the left half plane. This is indeed true in the large N limit of our 
matrix integral. 
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These identities are used along with the large N factorization 

(trA • trB) ~ (trA) (trB) (A.4) 
to derive various relations among disk correlators. 

A.l Loop equation for the resolvent 

An equation for the resolvent (4.11) is obtained if we take F = W(x). The identity (A.l) then gives 

1 n 1 

W\x) = - <tr[y'(X)W(x)] )-£- <tr[Y2W(x)] > . (A.5) 

' a=l P 

Then the identity (A. 3) applied to the last term gives 

1 I' dr 1 

<tr[Y2W(x)]) = - ^ — (dy a [W(x')W(x)Y a W(-x')]) 

= Lj) _(tr[W(x')W(x)]tr[W(-x')]) 
P Jim. 27ri 

JiR 2vrz x - x' 

Using the ^-product introduced in (5.6), the last line can be written as /3[W * The equation 

(A.5) can then be written as 

W{xf - V\x)W{x) + n[W * W](x) = ^ M ^M^ ) \ . (A.7) 

For a cubic potential the expectation value on the r.h.s. is a polynomial of degree one. Using an 
important property of the ^-product 

[/ * g] (x) + [g*f] (-x) = f(x)g(-x) , (A.8) 

which can be proved by deforming the contour of integration, one obtains a loop equation [31], which 
is a quadratic functional equation for W(x). The term linear in W can be eliminated by a shift 

, v def , 2V'(x)-nV'(-x) 

w(x) = W(x) — tt 1 -• (A.9) 

4 — n z 

The loop equation for w(x) is given by (5.3). 

A.2 Loop equations for the boundary two-point functions 

The boundary two-point functions of L-leg operators satisfy the recurrence equations 

D^ +l = W*D^\ D® +1 = W*D® (L>1) 

which can be derived by applying (A.3) with F = WS L+ iHS L . By applying (A.3) to F = WY„H 
and WHY a , we find 

= (XmD^ +W)*D , 
X — (tr[WHYf 1} ]) = D * (X (1) D^ + W), (A.10) 



fin 

32 



(2) 

and a similar pair of equations for Dq and D\ . The first equation of (A. 10) can be rewritten into an 
equation for the discontinuity along the branch cut, 

DiscL»{ 1} (x) = Z?o(-^)Disc(A(i ) L>S 1) (x) + W(x)), (A.ll) 
which implies that the recurrence equation can be extended to L = by defining 

n (l) def Dp n (2) def Dp 



Mi) Dp — 1 ' A(2)Do — 1 

Also, by applying (A. 1) to F = WH one finds 



(W + H)D = V'(x)D - P(x) - - (tr[WH £ Y 2 {a) ] ) , (A.12) 



where P(x) is defined in (5.12). 

By using (A. 8) to combine the two equations in (A. 10) and noticing that 



tr[WH ]T \ a) Yf a) ] = tr[H + (y-x)WH-W], 



a=l,2 

one finds a quadratic relation 

(y-x)D (x)+H -W(x)+ V)n (o) ^ a) ( 

a=l,2 
>(«)/ 



A>(z) ^ (A 2 a) nW^ Q) (-x) + A (a) n(")^(-x)) . (A.13) 



a=l,2 

By inserting the second of the equation (A. 10) into (A.12) one finds another quadratic equation, 

,(«) 

a=l,2 



(W(a;)+H)Do(ar)+-P(x)-'^ , (a:)-Do(z)- ^ nC«)£)p(-x) 

«=1,2 

= -D (s){ J] A (a ,n(-)Dj a) (-^) + nW(-x)}. (A.14) 



a=l,2 

These two equations are equivalent to (5.10) and (5.11). 

B 2D Liouville gravity 

In 2D Liouville gravity formalism, the c < 1 matter CFTs are coupled to the Liouville theory with 
central charge 26 — c and the reparametrization ghosts. For Liouville theory, we denote the standard 
coupling by b and the background charge by Q = b + 1/6. The central charge is given by 26 — c = 
1 + 6Q 2 . In our convention, b is always smaller than one. When the matter CFT is (p, q) minimal 
model, we have 

& = c=l-fc< (B.l) 

V i pi 

In studying the 0(n) model we used the parametrization n = 2 cos ir9. If 6 = l/p with p G Z, the 
model describes the flow between + 1) and (p— minimal models corresponding respectively 
to the dilute and dense phase critical points. 
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The primary operators in (p, q) minimal model fit in the Kac table which has (p — 1) rows and 
(q — 1) columns. The operator s has the conformal dimension 

(rg _ sp) 2_ (g _ p) 2 {r/b _ sb) 2_ (1/b _ b) 2 

K > s ~ 4p-q " 4 ' (B - 2) 

and are subject to the identification <j? rs = $ p _ r (? _ s . In 2D Liouville gravity coupled to the (p, q) 
minimal CFT, the operator $ r s is dressed by the Liouville exponential e 2otr ' s ^ or e ar - B< ^ depending 
on whether it is a bulk or boundary operator, so that the total conformal weight becomes one. This 
requires 

cx r ,s(Q - QV,s) + h r ,s = 1, 2av jS = Q ± (r/b - sb) . (B.3) 
B.l Conformal weight and scaling exponents of couplings 

In making the comparison between the matrix model and Liouville gravity, we start from the fact that 
the resolvent w(x) and its argument x correspond to the two boundary cosmological constants p B and 
p B . They couple respectively to the boundary cosmological operators and e^/ 6 , and are therefore 
proportional to p2 and p 1 ^ 1 , where p is the Liouville bulk cosmological constant. Since w(x) ~ x 1±e 
in the dilute and dense phases, we find 

(dilute): x = p B , w(x) ~ x 1+e = p B , b=(l + 6)- 1 / 2 , 

(dense): x = p B , w(x) ~ x 1 - = fi B , b = (1 - 6) +1 l 2 . 1 ' ' 

Suppose a boundary condition has a deformation parametrized by a coupling which scales like 
x p . Then the corresponding boundary operator should be dressed by the Liouville operator with 
momentum a = pb in the dilute phase and a = p/b in the dense phase. Then using (B.3) one can 
determine the conformal dimension of the operator responsible for the boundary deformation. Let us 
apply this idea to the deformations parametrized by p B , t B and A in Sect. 6. 

Deformation by /ig. p B scales like x l+e in the dilute phase and x 1 ^ 9 in the dense phase. So 
the corresponding operators are dressed by e^ 1+e ^ = e^' b in the dilute phase and e^ l ~ e ^/ b = e b ^ 
in the dense phase. The matter conformal dimension is zero, i.e., the operator responsible for the 
deformation is the identity $1,1. 

Deformation by t B scale like x 6 in the dilute phase and x~ 6 in the dense phase, so the cor- 
responding operator gets dressings with Liouville momentum b9 = \ — b in the dilute phase and 
— 8/b = b — I in the dense phase. We identify the operator with $13 (relevant) in the dilute phase 
and $31 (irrelevant) in the dense phase. 

Deformation by the anisotropy coupling A. The coupling scales like x l ~ e in the dilute phase. 
The corresponding operator is dressed by the Liouville momentum 6(1 — 9) = 2b — \ and therefore 
identified with $3 3. 



B.2 Conformal weight and scaling exponents of correlators 

After turning on the gravity, correlators no longer depend on the positions of the operators inserted 
because one has to integrate over the positions of those operators. The dimensions of the operators 
therefore should then be read off from the dependence of correlators on the cosmological constant p. 



34 



If we restrict to disk worldsheets, the amplitudes with n boundary operators Of dressed by Liouville 
operators and A bulk operators Oi dressed by e 2a ^ scale with p as 

(Of--- B O l ■ ■ ■ O m ) oc ^W-aE^-EA). (b.5) 

Suppose a disk two-point function of a boundary operator O b scales like x p . Then O b should be 
dressed with Liouville momentum (3 = \{Q — bp) in the dilute phase and (3 = |(Q — p/b) in the 
dense phase. We thus identify O b with the (r, s) operator in the Kac table if 

± p = r(l + 9) — s (dilute phase), ±p = r — s(l — 0) (dense phase). (B.6) 
B.3 Boundary two-point function 

In Liouville theory with FZZT boundaries, the disk two-point structure constant is 

(^ 7 (fe>^ )( Q- 2/3)/2bG(Q _ 2/3)G -l (2/? _ Q) 

a\p\z,s) s ^ + it + i8 ) S (p + it _ i8 -) S (p_ it + i8 ) 3 (p_ it _ ia y <■ ■ > 
Here the special functions G and S satisfy 

S(x + b) = 2sin(vr6x)S(x), G(x + b) = -^b^- bx T{bx)G(x), (B.8) 

V27T 



and similar equations with b replaced by 1/6. The boundary parameter s is related to the boundary 
cosmological constant by 

Pb = ( — ^To ) 5 cosh(2^6s), p B = (—^—) 1 cosh(27rs/6). (B.9) 
where p is the bulk cosmological constant and p is its dual, 

(p7T 1 (b 2 )b 2 ~ 2b2 ) 1 / b = {~p^{b~ 2 )b~ 2+2b ^) b . (B.10) 

Our loop equation (5.14) can be compared with 

d((3\t,s)d(Q -13- \\t ± f ,s) = F(/3){ cosh 2irb(t±if3) -cosh 2tt6s}, 



or the one with b replaced by 1/b. Similarly, the recurrence equation (5.7) among disk correlators of 
L-leg operators can be compared with 



d({3\t+%,s)-d{P\t-%,s) 
d(/3 + ±\t,s) 



GQ3){ cos27ri6(t + f ) - cos2vri6(t - § )}, 



rrm drf /~^~ r(i + 6 2 )T(1 - 26/?) 

G(/?) - "Vsinvr^ T(2 + b 2 -2bf3) ' (B>12) 
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C Properties of the function V r (r) 

The function V r (r) denned by 

, „ ( v def 1 P d W 
lo g y r (r) = - - / — 

^ ./-oo w 



e lWT sinh(7rru>) rc? 
sinh(-7ra;) sinh ^ ttuj 



(CI) 



satisfies the shift relations 



V r+1 (r) = 2cosh{ e{T± 2 i7Tr) ) 



K + i(t) = 2cosh(^^)K.(TTi), (C2) 



' t ± Z7rr 
"2 

which follow from the integration formula 

2uj 



log (2 cosh t) 



+ 



sinh(7ru;) ttuj 



(C.3) 



By deforming the contour of integration and applying the Cauchy theorem we can write the inte- 
gral (C.l) as the following formal series which makes sense for Re[r] > 0: 



lnK(r) = - f tT e ~nr _ y H! e ^^. (C4) 
2 ' n sm(7rn/y) ' n sm(mrU) 

n=l n=l 

The expansion for Re[r] < follows from the symmetry y r (r) = V t (—t). The expansion of the 
function V r at infinity is 



V r (r) 
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